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Spatial symmetries of the densities appearing in the nuclear Density Functional Theory 
are discussed. General forms of the local densities are derived by using methods of 
construction of isotropic tensor fields. The spherical and axial cases are considered. 

1. Introduction 

The nuclear Density Functional Theory (DFT) has become a conventional ap- 
proach to describing complex nucleiP^El In nuclear theory, the nuclear DFT is 
closely linked with the Hartree-Fock-Bogoliubov (HFB) approximation to the nu- 
clear many-body problem formulated in the position space. The mean value of the 
system's energy can here be expressed as a functional of the single-particle density 
matrices in the particle-hole (p-h) and particle-particle (p-p or pairing) channels. 
The density matrices in both channels are expressed through the scalar-isoscalar, 
scalar-isovector, vector-isoscalar and vector-isovector densities, which are nonlocal, 
i.e., they depend on two position vectors. In the local theory, e.g., the local den- 
sity approximation, the nonlocal densities can be represented in the energy density 
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through local densities and their derivatives. 

The nuclear self-consistent mean field resulting from the DFT may spon- 
taneously break the symmetry of the original many-body Hamiltonian. The 
symmetry-breaking phenomenon allows for capturing important physics effects and 
for including essential correlations in the many-body system. The existence of a 
self-consistent symmetry (SCS), i.e., a symmetry obeyed by the mean-field Hamil- 
tonian, reflects physical properties of the system. Due to the self-consistency, the 
SCS of the mean field is also the symmetry of the density matrix. This has a number 
of important consequences. First of all, to understand SCSs of DFT, it is sufficient 
to study symmetries of the density matrix without considering the specific form of 
the density functional and resulting mean-field Hamiltonian. We also realize that if 
the initial density matrix used at the first iteration of the self-consistent procedure 
has a certain SCS, then this symmetry will propagate through to the final DFT 
solution. This means that the choice of the initial density matrix is crucial for the 
proper physical content of the DFT results. Sometimes, e.g., to simplify numerical 
calculations, a SCS of the initial density matrix is ad hoc assumed. Such an as- 
sumption may lead to overestimation of the DFT energy and incorrect prediction 
of many properties of the system. Finally, assuming a SCS, one should know the 
general form of the density matrix possessing that symmetry. 

Some point symmetries and the associa ted s ymmetry-breaking schemes of the 
p-h densities have been investigated in Refs.EEl A detailed analysis of SCSs of the 
DFT densities in both channels has recently been presented in RefP, which will be 
referred to as I in the following. The present paper is a comment to I, with the main 
goal of detailing the consequences of spatial symmetries discussed therein. Although 
our discussion remains rigorous, in this contribution we simplify the notation so as 
to make the results more transparent and understandable. 

As in I, we treat the nonlocal and local densities as isotropic tensor fields, i.e., 
functions of the position vector (s) indep ende nt of other tensor quantities (material 
tensors) 3 Under this assumption, in RefsP^we introduced the Generalized Cayley- 
Hamilton (GCH) theorem for tensor fields, and we used it as a tool to study general 
forms of local densities. Since many different generalizations of the Cayley-Hamilton 
theorem exist in the literature, in Sec. [6] we give a brief explanation of the present 
version. 

In Sec. [2] we recall definitions of the density matrices and nonlocal and local 
densities. The spherical symmetry with and without space inversion is discussed in 
Sec.[3l In Sec.|4j we consider the axial symmetry with and without space reflection. 
The summary of our study is contained in Sec. El 

2. Density matrices and densities 

Within the HFB approach to the nuclear many-body problem, the mean value of a 
many-body Hamiltonian is a functional (the energy functional) of the p-h and p-p 
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density matrices denned, respectively, as 

p(rst,r's't') = (^lo^Orrtl*), (1) 
fast, r's't') = 4sV<#|<v_ s <_ t <a r . st |#), (2) 

where a^ st and a rs t create and annihilate, respectively, nucleons at point r, spin 
s=±i, and isospin t=±i, while \^f) is the HFB independent-quasiparticle state. 
The properties of density matrices that directly result from their definitions are the 
following: 

p*{r st, r's't') = p(r's't',rst), (3) 
fast, r's't') = -16ss'tt'p(r' - s' -t',r-s- t). (4) 

The dependence on the spin and isospin variables in the density matrices can be 
easily separated by expanding in the spin and isospin Pauli matrices & ss > and 
itf, (k = 1, 2, 3), respectively: 

fast, r's't') = l(po(r,r')8 ss < + s {r,r') -ow)tftt' 

+ z^^'P^r') + s k {r,r') ■ & as >)f*,, (5) 

k 

fast, r's't') = \(p (r,r')5 ss > +s (r,r') ■& ss ,)5 tt > 

+ |^(^p fc (r,r') + s fc (r-,r')-^)^, (6) 

k 

where k = 0, 1, 2, 3. The spin-isospin components of the p-h (p^, s^) and p-p (pk, 
§k) nonlocal densities are functions of two position vectors r and r' and have the 
following symmetry properties that result from Eqs. ([3]) and 

Pk{r,r') = p* k {r',r), 

s k (ry) = st(r',r), (7) 

for k — 0, 1, 2, 3, and: 

p k (r,r') = TPk(r',r), 

s k (r,r') = ±s k (r',r), (8) 

where the upper sign is for k = (isoscalars) and the lower for k = 1,2,3 (isovec- 
tors). 

In general, the p-h and p-p density matrices transform differently under the 
single-particle unitary transformations. However, it is proved in I that for the spatial 
transformations such as rotations, space-inversion, etc., the transformation rules 
for the p-p density matrix are the same as those for the p-h density matrix. These 
rules are obviously the same for all the isospin components. Therefore, in further 
discussion of the space symmetries, we shall omit the accent "breve" and the index 
k. Also, we shall not take into account conditions (O and jSJ, i.e., the hermiticity 
of the p-h densities and the symmetry or antisymmetry of the p-p densities. We 
mention only that the former condition ensures the reality of all the p-h local 
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densities, whereas the latter one results in vanishing of either the isoscalar or the 
isovector p-p local densities. 

Within the local density approximation, the energy functional is built from the 
local densities (r = r') and deriv atives thereof. The exact definitions of all used local 
densities are given, e.g., in RefsPES Here we only provide schematic definitions 
that clearly expose the corresponding spatial properties. The densities of interest 
are: 

• zero-order local densities 

— particle or pairing density 

p(r) = p(r,r') r=r , (9) 

— spin density 

s{r) = s{r,r r ) r = r , (10) 

• first-order local densities 

— current density 

j(r) = M(V-V>(r,r')]„ =r , (11) 

— spin-current density 

^^[(V-V')®^')]^,, (12) 

which is decomposed into the spin-divergence (trace of J) density J, 
the spin-curl (antisymmetric part of J) density J, and the traceless 
symmetric spin-current densitjQ J: 

•^) = M(V-V')- S (r,r')] r=r „ (13) 

J(r) = I[(V-V')x S (r,r')]_„ (14) 

J(r) = i [ (V - V) eg s(r,r' ) ] r=r , (15) 

• second-order local densities 

— kinetic density 

r(r)=[(V-V>(r,r / )] r=r , (16) 

— spin-kinetic density 

T(r)= [(V- V')s(r,r')] r=r , (17) 

— spin-tensor density 

F(r) = J[(V8V'+V'«V)..(r,r')] r=r ,. (18) 

We confine ourselves to the second-order derivatives as is usually done. But our 
analysis can also be extended to higher-order derivatives of the nonlocal densities 
that have recently been considered!^ 

a Underlined symbols stand for tensors after symmetrizing and subtracting the trace, e.g., J a() = 
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3. Spherical symmetry 

3.1. Nonlocal and local fields 

Let r and r 1 be two arbitrary, linearly independent position vectors. Then the vector 
product r x r' is the third linearly independent vector, and all three form a basis in 
the three-dimensional space of positions. The scalar products: r r — r 2 , r' -r' = r' 2 , 
and r -r' form three independent scalars quadratic in r, r'. It is impossible to form 
a cubic scalar because r ■ (r x r') = and r' ■ (r x r') = 0. Six possible outer 
products of the three vectors in question form the following second rank Cartesian 
tensors: 

- three quadratic tensors — r <g> r, r' ® r', and r ® r'; 

the first two tensors are symmetric and their traces are r 2 and r' 2 , respec- 
tively; 

the vector antisymmetric part of the third tensor is r x r' and its trace is 
r ■ r'. 

- two traceless cubic tensors — r (g> (r x r') and r' (g) (r x r') with the 
vector antisymmetric parts equal to r x (r x r') = (r • r')r — r 2 r' and 
r' x (r x r') = — (r ■ r')r' — r' 2 r, respectively. 

- one fourth-order tensor (r x r') ® (r x r') which can be expressed as a 
linear combination of the quadratic tensors with scalar coefficients: 



where (l) a h = <5 a 6 (a, b = x, y, z) is the unit tensor. 

Having listed all the independent scalars, vectors, and tensors that can be con- 
structed from vectors r and r', we are able to give general expressions for the 
nonlocal isotropic fields depending on the two position vectors. We note that (i) 
any scalar field must be an arbitrary function of the independent scalar functions: 



(ii) any vector field must be a linear combination of r, r', and r x r' with scalar 
coefficients: 



Q(r,r') = qn{r 2 ,r • r',/ 2 )r + q 12 {r 2 ,r ■ r' \r' 2 )r' + q 13 (r 2 ,r ■ r',r /2 )r x r'; (21) 



(r x r') ®{rx r') = (r ■ r')(r ® r' + r' ® r) - (r • r') 2 l 
- r 2 (r' ® r') - r' 2 (r ®r) + r 2 r' 2 l, 



(19) 



Q(r,r') = q (r 2 ,r-r',r' 2 y y 



(20) 
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and (iii) any symmetric traceless tensor field must be a linear combination of the 
five basic tensors: 

Q(r, r') = q 2 \ {r 2 , r ■ r' , r' 2 ) r <g> r + q 22 (r 2 , r ■ r', r' 2 ) r' <g> r' + q 2 3(r 2 ,r ■ r' , r' ) r ® r' 
+q2i{r 2 ,r ■ r', r' 2 )r ® (r x r') + q 25 (r 2 ,r ■ r' ,r l2 )r' ® (r x r') 
= q 2 i (r 2 , r • r',r /2 )(r ® r - \r 2 l) + q 2 2(r 2 ,r ■ r',r' 2 )(r' <g> r' - |r' 2 l) 
+g23(»- 2 ,f •r',r' 2 )[i(r ®r' + r'®r) - |(r-r')l] 
+g24(r 2 ,r • r',r' 2 )|(r ® (r x r') + (r x r') ®r) 

+g25(r 2 ,r • r' .r' 2 )\(r' <g> (r x r') + (r x r') ®r'). (22) 

In the expressions above, all g's are arbitrary scalar functions. Scalar fields always 
have the positive parity. The parities of vector and tensor fields are, in general, in- 
definite. However, since each independent vector or tensor field does have a definite 
parity, the vector and tensor fields of definite parities can be easily defined. 

It is readily seen from Eqs. (f20|) . ff2T|) . and (|22|) that the corresponding local 
fields, which depend on one position vector r = r' , only take very simple forms (cf. 
Appendix A in I): 

Q(r) = q (r 2 ), (23) 
Q{r)= qi (r 2 )r, (24) 
Q(r) =q 2 {r 2 ){r®r~\r 2 l). (25) 

3.2. Nonlocal and local densities 

3.2.1. Rotational symmetry SO(3) 

If we assume that the density matrices p and p, Eqs. ([T]) and @, are invariant under 
the three-dimensional rotations forming the SO (3) group, it immediately follows 
from Eqs. flS]) and ^ that the densities of type p are the SO(3) scalars while the 
densities of type s are the SO (3) vectors (note that the spin Pauli matrices are the 
SO(3) vectors). Therefore, the nonlocal density of type p takes a form of Eq. (|2U| : 

p(r,r') = Q (r 2 ,r-r',r' 2 ), (26) 

while the nonlocal density of type s has a form of Eq. (|2T|) : 

s(r,r') = gn{r 2 ,r ■ r',r' 2 )r + gi 2 (r 2 ,r ■ r',r' 2 )r' + gi3(r 2 ,r ■ r',r l2 )r x r', (27) 

where go, gn, g\ 2 , and q\$ are arbitrary scalar functions. All local differential 
densities ([TTj) - (fl"%|) can be calculated by differentiating Eqs. (|2"6")l and (|2T|) . like it 
was done in I. But to establish general forms of all local densities, it is sufficient to 
realize that all of them are local isotropic fields with definite transformation rules 
under SO(3) rotations. These rules can be deduced from the definitions (|9|)- ([T8|) . 
We see that the local densities p(r), r(r), and J(r) are scalar fields and all take the 
form of Eq. (f23|) . Similarly, the densities s(r), j(r), J(r), T(r), and F (r) are the 
SO(3) vectors; hence, are given by the form of Eq. (f24| . Finally, the spin-current 
density J(r) is the traceless symmetric tensor of the form ([25]) . 
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3.2.2. Rotational and mirror symmetry 0(3) 

When the density matrices are also invariant under mirror reflections, it follows 
from Eqs. ([5]), (JB]), and (|2"6"|) that p(r, r') should have positive parity. The spin Pauli 
matrices form an 0(3) pseudovector and thus the nonlocal spin density should be 
a pseudovector as well. On the right-hand side of Eq. (|2"T|) . only the last term is 
a pseudovector. Therefore, in the case of the 0(3) symmetry, the spin nonlocal 
density takes the form: 

s(r,r') = g((r 2 ,r -r',r' 2 )(r x r'), (28) 

meaning that the local spin density s(r) = 0. It is impossible to build a pseudovector 
from one vector r; therefore, local densities T(r) and F(r), being pseudovectors, 
must vanish, as well as pseudoscalar J(r) = and pseudotensor J(r) = 0. On the 
other hand, vectors j(r) and J(r) do not vanish and take the form ([24| . 

4. Axial symmetry 

Let us suppose that rotations and mirror rotations around one axis (say z-axis) 
are SCSs. This symmetry group will be denoted as O zX (2) cO(3). It is the direct 
product O zl (2) =Sz<g>SO J -(2) of the SO ± (2) group of rotations around the z-axis 
and the two-element group S z consisting of the reflection in the plane perpendicular 
to the symmetry axis and the identity. To investigate the axial symmetry, it is 
convenient to decompose the position vector r into the components parallel and 
perpendicular to the symmetry axis: 

r = r z +r ± , (29) 

which have different transformation properties under the O z_L (2) transformations. 
The component r± is a SO" L (2) vector whereas r z is not affected by the SO _L (2) 
rotations; hence, it is invariant under SO _L (2). On the other hand, r z changes its 
sign under the reflection S z , while r± is S z -invariant. 

4.1. Two-dimensional rotational symmetry SO 1 - (2) 

If the density matrices p and p possess the SO" 1- (2) symmetry, it follows from Eqs. ([5]) 
and ([6j) that all densities of type p should be SO x (2) scalars. Now it is convenient 
to decompose densities of type s into components perpendicular and parallel to the 
symmetry axis: 

s(r,r') = s ± (r,r') + s z (r,r'). (30) 

The perpendicular components s±(r,r') are SO _L (2) vectors and the parallel com- 
ponents s z (r,r') are SO _L (2) invariants. 

The following SO -1 " (2) scalars can be constructed from the perpendicular and 
parallel components of the position vectors r and r': z, z', the z-coordinates of r z 
and r' z , respectively, r\ = r± ■ r±, r± ■ r' ± , r' 2 = r' ± ■ r'j_, and r z ■ (r_i_ x r' ± ) and 
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r' z ■ (r± x r' x ) = (z' / z)r z ■ (r x x r' x ). As we see the two latter scalars are dependent 
on each other because r z \\ r' z . The square (r z ■ (r x xr' ± )) 2 = z 2 (r x r' 2 — (r x -r' ± ) 2 ) is 
dependent on other scalars, namely on z 2 , r x , r' 2 and r±-r' x . Hence, the nonlocal 
scalar densities have the form 

p(r, r') = g (z, z', r 2 x ,r x -r' x ,r'l,r z ■ (r± x r' x )), (31) 

where go is an arbitrary scalar function, linear in the last argument. There are three 
vectors invariant under SO x (2) and parallel to the symmetry axis, namely r z , r' z , 
and r x x r' x . Consequently, the z-component of s(r, r') is of the form 

s z (r,r') = g z (z,z',r 2 _,r± ■r' x ,r'l,r z ■ {r± x r' ± ))r z 
+ Q' z ( z > z '> r l> r -L -r' ± ,r'l,r z ■ (r x x r' ± ))r' z 
+ <p z (z,z',r ± ,r x -r' ± ,r'l,r z ■ (r x x r' x ))(r x x r' x ), (32) 

where g z , g' z , and ip g are scalar functions. There exist four independent SO J_ (2) 
vectors that lie in the plane perpendicular to the symmetry axis: t*Xi v\, r z x r± 
and r' z x r' x . Hence, the SO J -(2) vector component of the nonlocal spin density 
takes the form: 

s±(r,r') = g±(z,z',r ± ,r± -r' x ,r'l,r z ■ (r± x r' ± ))r± 
+g' x {z,z',r x ,r ± -r' x ,r'l,r z ■ (r x x r' x ))r' x 
+ip x (z,z',r x ,r x -r' ± ,r x ,r z ■ (r x x r' ± ))(r z x r x ) 
+^ ± (z,z',r 2 x ,r x -r' ± ,r' 2 ,r z -(r x xr' ± ))(r' z xr' x ), (33) 

with g x , g' x , ip±, and ip' x being scalar functions. General forms of the local zero- 
order densities are obtained from Eqs. ([3Tj) . (f32|) . and (|33|) by putting r z — r' z and 
r x = r' x . This gives: 

P (r) =p (z,r x ), (34) 
s z (r) = p z (z,r x )r z , (35) 
s x (r) = p±(z,r x )r x + (/>±(z,r x )(r z x r x ), (36) 

where po, p z , p±, and 4> x are arbitrary scalar functions of z and r x = v/rT. The 
general form of the remaining scalar densities r(r) and J(r) is the same as that of 
Eq. {3U). 

The components of the differential operator, 

V = V z + V x , (37) 

have the same transformation properties under SO -1 " (2) rotations as the correspond- 
ing components of the position vector (|29|) . Hence, the densities j(r), J(r), T(r), 
and F(r) have the components parallel to the z-axis that are invariant under 
SO x (2) and the SO x (2) vector components that are perpendicular to the sym- 
metry axis. Therefore, they all take general forms given by Eqs. (|35|) and (|36p . 

It follows from the definitions (fT2|) . (fl5|) . and (|37|) that the components {S} az 
(a = x, y) of the symmetric traceless spin-current density form the SO x (2) vector 
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while the components (J) a & (a, b = x, y) form the 80^(2) symmetric traceless 
tensor. The following four symmetric traceless tensors can be formed with vectors 
r z and r±: 

l{r z ®r±+r x ®r z ), (38) 



2 



r z ® (r z x r ± ) = \(r z <g> (r z x r ± ) + (r z x r ± ) ® r z ), (39) 
rj_ ® rj_ = r_L 8) rj_ - |rjlx, (40) 
rj_ (r z x = ±(rj_ ® (r z x r±) + (r z x r±) <g> rj_), (41) 

where (lj_) a & = (1 — S az )S a b (a, b — x, y, z) is the unit tensor in the plane 
perpendicular to the symmetry axis. Two tensors (|38[) and (1391) transform under 
SO -1 - (2) like vectors perpendicular to the z-axis. The tensors (|40]) and (|41[) are the 
SO -1 - (2) tensors. Consequently, the general form of the symmetric traceless spin- 
current density is: 

J(r) = z ±(z, r±) r z ® r± + 6 zz ±(z 7 r±) r z g (r z x n) 

+ fl_u_(.z,r_i_) r_i_ igiTx + ^± z ±(z,rj_) rj_ g (r z x rjj , (42) 

where 6* z j_, # zz j_, #x± ; an d #_l z _l are scalar functions. We note in passing that in 
the case of the axial symmetry, all local densities formally look like the spherical 
symmetric nonlocal fields (j2"0)) . (|2ip . and (|2"2"|) with r = r z and r' = r±, provided 
that r z ■ r± = 0. 



4.2. Two-dimensional rotational and mirror symmetry O z± (2) 

It follows from the S z invariance of the density matrices Q and © that p(r, r') 
and s z (r, r') do not change signs under S z , while s±(r, r') does. All scalar functions 
invariant under S z depend on z, z' only through z 2 — r z ■ r z , z' 2 — r' z ■ r z , and 
zz' = r z ■ r' z . For instance, Eq. (|3"Tj) for the nonlocal scalar density now takes the 
form 

p(r,r') = p (2 2 ,zz',z' 2 ,ri,r ± -r'^rf). (43) 

There exist two pseudoscalars formed from r z , r' z , r±, r' ± , namely r z ■ (rj_ x r' ± ) 
and r' z ■ (r± x r' ± ) = (zz' / z 2 )r z ■ (rj_ x r' ± ) equal to each other up to the scalar 
factor zz' jz 2 . The square (r z ■ (r± x r ± )) 2 = z 2 (r\r' 2 — (r± ■ r'j_) 2 ) is, of course, a 
scalar. 

To fulfill the transformation rules, the general forms of Eqs. ([3^)) and (f3"3")l of the 
components of the nonlocal spin density should be modified in the following way: 

s z (r,r') = (r z ■ (r x x r' ± ))g z (z 2 , zz', z' 2 , r\, r ± ■r'^,r'l)r z 
+ (r' z ■ {r x x r' ± ))g' z (z 2 ,zz',z' 2 ,r 2 ± ,r ± -r' ± ,r'l)r' z 
+<p z (z 2 ,zz',z' 2 ,r 2 ± ,r±-r' ± ,r'l)(r ± xr' ± ), (44) 
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and 

s±(r,r') = (r z ■ (r± x r' ± ))g ± (z 2 1 zz', z' 2 , r ± , r± ■ r' x ,r'l)r± 
+(r' z ■ (r x x r' ± ))g' ± (z 2 , zz' , z' 2 , r\, r x ■ r' x> r'l)r' x 
+L P± (z 2 ,zz',z /2 ,r 2 L ,r ± ■r' x ,r'l){r z x r x ) 
+^ ± (z 2 ,zz',z' 2 ,ri,rx-r> ± ,r' 2 )(r' z xr' ± ). (45) 

It follows from Eqs. (|43p . (|4"4")) . and (|4"5")) that the local zero-order densities for 
r z = r' z and r x = r' ± can be written in the general form: 

p(r) = Po (z 2 ,rl), (46) 
s z (r) = 0, (47) 
s ± (r) - q> ± {z 2 ,ri)(r z x r ± ). (48) 

The local kinetic density r(r) is of the form (|4"6")l too. On the other hand, the pseu- 
doscalar density J(r) vanishes. The densities T{r) and F(r), are pseudovectors; 
hence, they take the form (|48|) . On the other hand, vectors j(r) and J r (T > ) are linear 
combinations of the components of the position vector: 

j(r) = L z (z 2 ,rl)r z + ix(z 2 ,rl)rx. (49) 

The spin-curl J(r) takes a similar form to that of (|49)). Finally, J(r) is a pseudoten- 
sor. Therefore, as follows from (|25[) . its general form is given by: 

i(r) = 6* zz _L(z 2 ,ri) r z g (r z x r ± ) 

+ Oxzx (z 2 ,r 2 ± ) r ± <S> (r z x r_Q . (50) 

5. Summary 

In the DFT, for both theoretical and practical reasons, it is important to know 
what general forms of densities are which obey SCS of interest. In the case of the 
space symmetries, such general forms can be established by means of methods of 
constructing the isotropic tensor fields. 

For the spherical symmetry, the local densities are the isotropic scalar, vector, 
or (the second rank) tensor fields, depending on the position vector r. The form of 
an isotropic field with given rank is unique and determined through one arbitrary 
scalar function. In particular, the parity of the field is unique for a given rank. 
Pseudoscalar, pseudovector, and pseudotensor fields do not exist. This is why in 
the case of the rotational and mirror symmetry, the pseudoscalar (spin-divergence), 
pseudovector (spin, spin-kinetic and tensor-kinetic), and pseudotensor (symmetric 
spin-current) local densities vanish. 

For the axial symmetry, the local densities are isotropic fields depending on two 
components of the position vector: r z and rx- The case of SO -1 " (2) is interesting as 
it allows us to better understand the assumption of the isotropy of a field. Indeed, in 
this case it might seem that the local densities are fields, depending on the 80^(2) 
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vector r±. However, such fields are not isotropic. There is another vector, r z , which 
plays the role of a material vector fixed by the direction of the symmetry axis of 
the system. 

6. Appendix: The Generalized Cayley-Hamilton Theorem 

The original Cayley-Hamilton The orem s tates that every square matrix satisfies its 
own characteristic equation (cf. e.gJ ^^^ J). It immediately follows from the theorem 
that the second rank SO (3) Cartesian tensor q satisfies the equation 

- q 3 + giq 2 - g 2 q + <73 = 0, (51) 

where q\ is the trace, qi is the sum of the principal subdeterminants, and is the 
determinant of q. Hence, a second rank tensor field Q(q) being the power series of 
the tensor q 

Q(q) = col + c x q + c 2 q 2 + c 3 q 3 + . . . (52) 

can be sumed up to the form: 

Q(q) = Po(<7i,<72,<73)l + pi(<?i,g 2 ,<73)q + p 2 (gi, g 2 , <j 3 )q 2 , (53) 

where p±, pi and ps are functions of scalars q\, qi and q%. The power series (|52|) 
is an isotropic function of q because it does not contain any other tensorM- 1 For a 
symmetric traceless (qi = 0) tensor Q (a quadrupole tensor) Eq. (|S*3"j) takes the 
form: 

Q(q) = £>i(<7 2 ,g 3 )q + £2(92,93)9^ (54) 

with two functions Qi and g 2 of the two scalars g 2 and q%. Eq. (|54[) is a direct 
consequence of the Cayley-Hamilton Theorem and thus can be called the Cayley- 
Hamilton theorem for the quadrupole tensors. It can be generalized for the isotropic 
tensor field (q^ , q<- y ^ , . . . ) with an arbitrary multipolarity L being a function 
of one or a few spherical tensors q^ x \ q^ x \ ... of ranks A, A', . . . , respectively. 
The Generalized Cayley-Hamilton Theorem has the following general form: 

QW( g W,^^..0 = fc( ^'''' ) J R,(g)TW( g W,^'),...), (55) 

k=l 

where q stands for the set of independent scalars, are some definite fundamental 
tensors, all constructed from q( x \ q( x \ . . . , and Rk are arbitrary scalar functions. 
The number k(L, A, A', . . . ) depends on the ranks of the all involved tensors. To 
find the number and the forms of the fundamental tensors in a general case may 
appear to be a difficult task. A systematic method of constructing them in the case 
of the tensor fields depending on one tensor are presented in Refs!^^ Eqs. (f20| - 
(f22|) are examples of the GCH theorem for A = A' = 1 and L = 0, 1, 2, whereas 
Eqs. (f23| -(|25 |) — for A = 1 and L — 0, 1, 2. In t he nuclear collective model the 
GCH theorem was used for A = 2 and A = 3.EHH 
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